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ABSTRACT 



Values of system availability are computed for a system 
whose failure density is exponential (X) and whose repair 
density is special Erlang (y,a). The system is modeled as 
an alternating renewal process. The Laplace transform of 
the availability function is developed and inverted (by 
using a numerical procedure) to obtain the values of the 
system availability. 

Corresponding values of the long-run average system 
availability and the average availability over the first 
mission time are also computed. Comparisons are made 
which establish that using the long-run average value to 
approximate the true availability over the first mission 
time is a conservative practice. 
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I . INTRODUCTION 



A . THE SYSTEM 

Most physical systems people encounter are of the type 
which operate for a time, fail, are repaired and put back 
into operation. The system then may be said to have two 
states. Truely, each of these states may be the union of 
numerous sub-states, but the two major states prevail. 

Some usual terms applied to the first state are 
'operational', 'available' or 'up'. This is the state or 
condition of the system when it is capable of accomplishing 
the task for which it was designed and built. The second 
state, the compliment of the first, is the less desirable 
state and is usually termed 'failed', ' nonoperatlonal ' , 
'unavailable' or 'down'. 

Hereinafter the system states will be referred to as 
' up ' or ' down ' . 

B. DEFINITION OP AVAILABILITY 

A good deal of effort is being expended by the U.S. 
Governmient and by private Industry to obtain measures of 
certain quality characteristics applicable to such systems 
as is described in A. Hosford [2] has discussed several 
of these characteristics. 

One such quality characteristic that is of primary 
concern is system availability. The 'system availability' 
or simply 'availability' is the probability that the 
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system is up at any time t^O, given it was up at time 
zero . 

The availability function, n(t), is a real valued 
function which gives this probability for any time t ^0. 

C. METHOD CURRENTLY PRACTICED TO DETERMINE AVAILABILITY 
In reliability literature, availability is often 
defined differently as: 



MTTP + MRT 

where MTTP is Mean-Tlme-To-Pallure and MRT is Mean-Repair- 
Time. This definition provides the only method currently 
in wide spread use for determining availability. 

Note: This quantity has the variable name ZBAR in the 

computer program and output (Appendices A and C). 

Cox [1] has shown that regardless of the distributions 
of uptimes and downtimes 

11m Ti(t) = A. 
t oo 

That is, A is actually the steady-state system availability. 

Similarly, it may be shown that the long-run average 
availability 

1 ^ 

11m -r f n ( x)dx = A . 
t^«> 0 
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Jacobson [ 3 ] describes a statistical procedure for 
obtaining a confidence interval estimate for A. 

D. OBJECTIVES OP THIS RESEARCH 

The objectives of this research were the following: 

1. To compute exact values of the availability function 
for a system with given distributions of uptimes and 
downtimes . 

2. To compare true availability of a system with that 
obtained by methods currently in practice. 
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II. SUMMARY AND CONCLUSIONS 



A physical system having an 'up' state and a 'down' 
state and with the ability to transit between the two 
states was modeled as an alternating renewal process for 
the purpose of finding the system's availability function, 
n(t). An explicit general expression for n(t) was not 
found because of the complex structure of the model 
equations under the assumed distributions of uptimes and 
downtimes. However, n (s), the Laplace transform of n(t), 
was developed from the assumed probability density functions 
and the renewal density functions of the model. 

A computer program was written to use a numerical 
procedure to Invert the Laplace transform for specific 
values of the parameters and of time, the Independent 
variable . 

Upon consideration of the results and some sample 
problems (see Appendices A and B) the following conclusions 
were made. 

1. Under the fairly realistic assumptions of the 
model, exact values of n(t) can be and were computed. 

2. Use of long-run average availability to approximate 
the true availability over the first mission time, as is 
currently being done in practice, is a conservative 
procedure . 
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3. The amount by which the long-run average availability 



underestimates the true availability during the first 
mission time is dependent upon the long-run value, i.e., 
the lower the long-run value, the more it underestimates 
the true availability. 
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III. THE MODEL 



A. DESCRIPTION 

The system was modeled as an alternating renewal process 
with two states. That is, the new system was assumed to 
start operating in the up state. The time spent in this 
state before the first failure, whether in use or just 
available for use, was a random variable U^. Upon failure 
the system entered the down state. The time spent down 
before the first repair, whether awaiting repairs or 
actually being repaired, was a random variable D^. The 
length of the first cycle was the random time + D^ . 

It was assumed the system kept alternating between up 
and down and each repair returned it to 'like new' condi- 
tion so that the uptimes U^, 1=1,2,..., were taken to be 
Independent, identically distributed random variables. 

They were assumed to be exponentially distributed with 
parameter A >0 and probability density function 



Note: All time units in this thesis are mission times. 

The downtimes D^,i=l,2 , . . . , were also assumed to be 
Independent, identically distributed random variables and 
independent of the uptimes. They were assumed to have a 




t < 0 



0 < t 
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special Erlang distribution with parameters y >0 and 
a=l,2,..., and with probability density function. 



rpCt) 



0 



y 

(a-D! 



, a-1 -yt 
t e 



t < 0 
0 < t 




The Special Erlang Density Function 



Weiss [5], through a simple argument, has shown the 
renewal density functions for such a renewal process to be 

t 

WpCt) = fyCt) + / W^(t) fy(t-x) dx 



t 

W^(t) = / Wq(x) rp(t-x) dx 
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where 0 and 1 are subscripts for ’up’ and ’down’ respec- 
tively and Wj(t) Is the renewal density for the event ’’end 
of stay in state J ” . Weiss has also given the availability 
function as: 

t 

n(t) = R(t) + / W,(t) R(t-x) dr 

0 ^ 

where R(t) is the reliability function, l.e., 

R(t) = / 

t ^ 

B. SOLUTION 

No straight forward solution of these equations for 
ri(t) in explicit closed form has been found except in the 
special case where fy and r^ are both exponential density 
functions. However, by finding the Laplace transforms of 
the probability density functions and the renewal density 
functions the Laplace transform of the availability function 
was found. 

1 . Laplace Transforms 

The Laplace transform of a continuous or sectlonally 
continuous function g(t), t in [0,<»), is 

* “ _o +- 

g (s) = / e g(t)dt ; s is a complex constant 
0 

whose real part is large enough that the integral converges. 
The Laplace transform of fy(t) is 
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f,j(s) = / e f.,(t)dt = / e A e dt 
0 0 



^ /-°° / 1 j. \ -(A+s)t A 

f (A+s) e dt = 



A+s 



A+s ' 



Similarly, rj^(s) may be shown to equal 



y+s 



a 



The Laplace transforms of the renewal densities are 



w*(s) = fyCs) + w*(s) fyCs) 



*/ \ */ \ * / \ 
w^(s) = Wq(s) rj^(s) . 



These Imply 



w!(s) = 

(A+s)(p+s)^ - 



Noting that In this case 



n(t) = 



WQ(t) 



the Laplace transform of the availability function was 
found to be 

„‘(s) = . 

( A+S ) ( y+s ) ^ - Ay“ 

2 . The Inverse Transform 

The desired solution for the availability function 
Is the Inverse transform of n (s). As before, no method 
was found to obtain an explicit general expression for n(t). 
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Therefore a numerical Inversion procedure was used to 
compute values of ri(t) for particular values of the 
parameters 1, y, and a and for values of t=.05, .10, 

.15, ...j 1.0 mission times. 

Several values were chosen for the parameters X, 
y and a over a range that was judged to be of practical 
Interest. Those values and the corresponding values of 
n(t) are shown in Appendix A. 



IV. DESCRIPTION OF THE COMPUTER PROGRAM 



The computations were done by a computer program written 
In FORTRAN IV programming language and run on the IBM 
system/360, model 67 computer at the Naval Postgraduate 
School (see program listing. Appendix C). The program 
listing contains enough explanatory comments to enable the 
reader familiar with FORTRAN to follow the procedure. 

The following are the major steps: 

1. The polynomial In the denominator of n (s) was 

expanded to the form aQ + a^s + a2S^ + ... + . That Is, 

(Ws)(yts)“ - = ? (“) (X+ . 

j = l ^ ^ 

Note that ag equals zero. 

2. A computer library subroutine which employs the 
Newton-Raphson successive approximations technique was used 
to find the roots of this polynomial. 

3. After the roots were checked to see that they were 
all different the following formula was used In computing 
the values of ri(t). Often called Heaviside's expansion 
formula [4] It Is an application of the residue theorem 
from theory of complex variables. 

Heaviside's Expansion Formula: let P(s) and Q(s) be 

polynomials where P(s) has degree less than that of Q(s). 

If Q(s) has distinct roots A^^, k=l,2 ,...,n, then the 

P ( s ) 

Inverse Laplace transform of is 
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L 



-1 




n 

Z 

k=l 



P(A^) A^t 



By application of this formula 



L“^(n*(s) ) 



n(t ) 



a+1 

k=l y+Aj^+a(A+A^) ^ 
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APPENDIX A 



RESULTS 



This appendix contains tabulated values of the avail- 
ability function for the indicated sets of parameters and 
the indicated time in mission-time units. 

The values of PZERO given are such that PZERO equals 
the probability of a repair in less than one mission 
time , 1 . e . , 



Z1 is the average availability over the first mission 
time . 




n(0)+n(l) . 
2 



19 

2 n(.05i)) 

1=1 



ZBAR is the asymptotic value of n(t) for large t. ZBAR 
is also the long-run average availability. 



ZBAR 



MTTF ^ 1/A ^ y 

MTTP + MTR ~ 1/A + a/y y + aA 



Example: For a system whose uptimes have the expo- 

nential (A = .005) distribution and whose downtimes have 
the Erlang (y = 2.03, a =2) distribution, the true system 
availability at a half mission time is n(.5) = .99777. 
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The true average availability over the first mission 
time is Z1 = .99792 and the long-run average availability 
is ZBAR = .99510. 
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ALPHA= 2.0 P7ER0 



TIME 


.005 


.010 


.050 


0.05 


.99975 


.99950 


.99751 


0.10 


.99950 


.99901 


. 99504 


0.15 


.99926 


. 99852 


.99263 


O 

• 

O 


.99902 


.99805 


.99027 


0.25 


.99879 


. 99759 


. 93799 


0.50 


.99357 


.99714 


.98580 


0.35 


.99836 


. 99671 


. 98369 


o 

• 

o 


. 99815 


. 99631 


.98167 


0.45 


.99795 


.99591 


.97975 


0. 50 


.99777 


. 99554 


. 97793 


0.55 


.99759 


.99519 


.97620 


0.60 


.99742 


. 99485 


.97456 


0.65 


.99726 


. 99454 


.97302 


0.70 


.99711 


.99424 


.07157 


0.75 


.99697 


. 99396 


. 97021 


0.80 


.99684 


. 99369 


. 96804 


0.85 


.99671 


.99344 


.96774 


0. 90 


. 99660 


.99321 


. 9666? 


0.95 


.99649 


.99299 


.96559 


1.00 


.99639 


. 99279 


.96462 


Z1 


.99792 


. 99584 


, 97945 


Z0AR 


.90510 


. 99024 


.95305 



0.60 


MIJ = 


2.03 




LAMBDA 


. 100 


. 200 


.300 


. 400 


99502 


.09007 


. 98514 


. 98023 


99011 


.93032 


.97063 


. 96103 


98531 


.97084 


.95658 


.94254 


08064 


. 96167 


.94307 


, 02484 


97614 


, 95286 


.93016 


.90802 


07130 


.94444 


.01788 


, 89210 


96766 


. 93642 


, 90626 


. 37712 


96370 


.92382 


.89530 


. 86308 


95995 


. 02164 


. 88501 


. 34997 


05639 


.91489 


.87533 


. 33778 


95303 


.90854 


.86640 


,82647 


94987 


. 90260 


. 85804 


.81602 


94689 


.89706 


.85029 


. 80639 


94410 


. 89189 


. 84312 


.79754 


94149 


. 88709 


.83650 


,78943 


93905 


. 88264 


.83040 


.78201 


93678 


. 87851 


. 82470 


.77524 


03466 


.87470 


. 81965 


.76908 


93270 


. 87119 


.31494 


.76349 


93087 


. 86795 


. 81064 


.75842 


95954 


.92151 


.88574 


.85203 


91031 


.83539 


.77136 


. 71731 
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ALPHA= ?.0 



PZERO= 0.70 
LAMBDA 



MU= ?.45 



TIME 


.005 


.010 


.050 


. 100 


.200 


.300 


.400 


0.05 


.99975 


.99950 


.99751 


.99502 


.99007 


.98515 


,98025 


0.10 


.99950 


.99901 


.99506 


. 99014 


. 98037 


.97071 


.96114 


0.15 


.99926 


.99853 


.99267 


.98539 


.97100 


.95682 


.94286 


0.20 


.99903 


.99806 


.99036 


.98082 


.96202 


.94360 


.92553 


0.25 


.99881 


.99762 


.98815 


.97646 


. 95349 


. 93109 


.90925 


0. 30 


, 99860 


. 99719 


.98605 


.97231 


.94544 


.91935 


. 89404 


0.35 


.09839 


.99679 


.98406 


.96840 


.03788 


. 90840 


, 87992 


C.40 


.99820 


.99641 


.98219 


.96472 


.93081 


.89822 


,86689 


0. 45 


.99802 


. 99605 


.98043 


.96128 


.92424 


.88881 


.85492 


0.50 


.99785 


.99571 


.97878 


.95807 


, 91816 


. 88016 


. 84397 


0. 55 


.99770 


.99540 


.97725 


. 95510 


.91254 


.87222 


.83400 


0.60 


.99755 


.99511 


.97582 


.95234 


.90738 


. 86497 


. 82405 


0.65 


.99741 


.99484 


.97450 


. 94979 


,90264 


. 85836 


. 81676 


0. 70 


.99729 


.99458 


.97328 


.94745 


.89831 


. 85236 


. 80938 


0.75 


.99717 


.99435 


.97216 


.94529 


.89436 


. 84693 


.80275 


O.PO 


.99706 


. 99414 


.97112 


.94331 


.89076 


.84203 


.79681 


0.85 


.99696 


. 99394 


.97017 


,94151 


. 88750 


. 83761 


.79150 


0.90 


.99687 


.99376 


.96930 


.93985 


, 88454 


. 83364 


, 78678 


0.95 


.99679 


. 99359 


.96850 


.93835 


.88186 


, 83008 


.78257 


1.00 


.99671 


.99343 


.96 777 


.93698 


. 87945 


. 82689 


,77885 


Z1 


.99803, 


.99606 


.98056 


.96170 


,92566 


. 89170 


. 85968 


ZBAR 


.99593 


. 99190 


.96078 


.92453 


.85965 


. 80328 


. 75385 



20 





ALPHA= 2.0 


PZERO 


1= 0.80 


MU = 


3,00 












LAMBDA 








TIME 


. 005 


.010 


.050 


. 100 


.200 


.300 


.400 


0.05 


.99975 


. 99950 


.99751 


. 99503 


.99008 


.90516 


,98027 


0.10 


.99951 


.09901 


.99508 


.99018 


.98045 


.97083 


.96130 


0.15 


.90927 


. 99854 


.99273 


.98552 


.97125 


. 95719 


,94334 


0. 20 


. 99905 


. 99809 


.99050 


,98109 


.96255 


.94437 


.92656 


0.25 


.90883 


.99767 


.98839 


.97692 


.95441 


. 93245 


. 01103 


0. 30 


.99863 


. 99727 


. 98642 


.97304 


, 94607 


.92146 


.00680 


0. 35 


.99845 


. 09690 


.98459 


.96944 


.93992 


.91141 


.88386 


0.40 


.99027 


.99655 


.98290 


.96613 


. 93357 


. 90226 


. 87216 


0.45 


.99812 


. 99624 


.98135 


. 96310 


, 92779 


.89400 


.86166 


0.50 


.99797 


.99595 


.97993 


.96034 


.92256 


. 88657 


.85227 


0. 55 


.90784 


. 99568 


.97863 


.95703 


.91784 


, 87991 


. 84393 


0. 60 


.99772 


.99544 


.97746 


.95556 


.91360 


.87398 


. 83655 


0.65 


.99761 


. 99522 


.97640 


.95352 


.90981 


.86871 


. 83004 


0.70 


.99751 


.99502 


.97544 


.95168 


.90643 


.06404 


.32433 


0.75 


.99742 


.9948^ 


.97450 


.95003 


.90342 


. 059Q3 


.01933 


o 

CD 

• 

O 


.90734 


.99468 


.97381 


.94856 


. 90076 


.85632 


. 81498 


0. R5 


. 99726 


. 99454 


.07311 


.94724 


.89840 


.8531 5 


.81120 


0.90 


.99720 


.994-41 


.97249 


.94607 


. 8963? 


,85037 


. 8070? 


0.05 


.99714 


. 99429 


.97194 


.94503 


. 89449 


. 84796 


. 80509 


1.00 


.99709 


. 99419 


.97145 


.94411 


.89288 


.84586 


. 80266 


Z1 


.09817 


.99635 


.98195 


.96442 


.93085 


. 89915 


. 86920 


ZBAR 


.99668 


.99338 


.96774 


.93750 


. 88235 


, 83333 


.78947 
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ALPHA= 2.0 PZFR0= O.^O MU= 3.90 



LAMBDA 



TIME 


.005 


.010 


.050 


. 100 


.200 


. 300 


, 400 


0.05 


,99975 


.99950 


.99752 


.99504 


.99011 


.98520 


.98031 


0. 10 


.99951 


.99902 


,99512 


.99026 


. 98061 


.97106 


.96161 


0.15 


.99928 


.99857 


.99285 


.98575 


.97171 


,95788 


. 94426 


0.20 


.99907 


. 99814 


. 99074 


. 98158 


.96352 


. 94582 


.92846 


0.25 


.99888 


. 99775 


.98881 


.97777 


,95608 


.93492 


,91428 


0.30 


.99870 


.99740 


.98707 


.97433 


, 94 <340 


. 92520 


.90170 


0.35 


.99854 


. 99708 


.98550 


.97125 


,94346 


.91661 


.89066 


0.40 


. 99840 


. 99680 


.98410 


,96851 


. 93823 


.90909 


. 88106 


0.45 


,<=9827 


. 99654 


.98287 


. 96610 


. 93364 


. 90255 


. 87276 


0. 50 


.99816 


. 99632 


.98178 


.96399 


.92965 


. 89689 


. 86565 


0.55 


.99806 


.99613 


.98083 


.96215 


.92619 


. 89204 


. 85958 


0.60 


.99797 


. 99505 


,98000 


.96054 


.92321 


. 88780 


.85444 


0.65 


.99790 


. 99581 


.97928 


.95916 


.92066 


.88435 


. 85011 


0.70 


.99783 


. 99568 


.97865 


.95796 


. 91847 


. 88136 


.84647 


0. 75 


,99778 


. 99556 


.97811 


.95694 


.91661 


.87884 


.84343 


0.80 


.99773 


.99547 


.97765 


.95605 


.91503 


. 87672 


. 84090 


0.85 


.99769 


.99538 


.97725 


.95530 


. 91370 


, 87494 


.83881 


0. 90 


.99765 


.99531 


.97691 


,95466 


.91257 


.87346 


. 83708 


0.95 


.99762 


.99525 


.97662 


.95412 


. 91162 


. 87223 


. 83566 


1.00 


.99759 


. 99520 


.97637 


.95365 


.91083 


.87121 


.83450 


Z1 


.99838 


. 99676 


.98399 


.96841 


.93850 


.91013 


. 88323 


ZPAR 


.99744 


. 99490 


.97500 


.95122 


.90698 


. 86667 


, 82979 
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ALI 


PHA= 2.0 


PZERQ 


= 0.95 


MU= 


4.80 












LAMBDA 








TIME 


.005 


.010 


.050 


.100 


.200 


.300 


.400 


0,05 


.99975 


.99950 


.9975? 


.99505 


.99013 


.98524 


. 98037 


0,10 


.9<3952 


.99903 


.99516 


.99035 


. 98080 


.97134 


.96198 


0.15 


.00930 


. 09859 


. 99298 


.08602 


.97225 


.95868 


.94531 


0.20 


.99010 


.99820 


.99102 


.98212 


. 96460 


.94741 


. 93057 


C.25 


.90892 


. 09784 


.98927 


. 97868 


.95787 


.93757 


.91776 


0.30 


.99877 


.99753 


.98775 


.97567 


.95204 


.92910 


. 90631 


0. 35 


.99863 


. 99727 


.98643 


.97308 


. 94705 


.92189 


. 89755 


0. ^0 


. 99852 


. 99704 


.98529 


.97086 


. 94282 


.91581 


. 88981 


0.45 


.00842 


.99684 


.98433 


.96899 


.93925 


.91074 


.88340 


0.50 


.09833 


. 09667 


. 98352 


. 96741 


.93628 


. 90654 


.87813 


0.55 


.99826 


. 99653 


.98283 


.96603 


.93381 


.90308 


.87333 


0.60 


.99820 


.99641 


.98226 


.96498 


. 93177 


. 90025 


. 87035 


0.65 


. 99815 


. 99631 


.98178 


.96407 


.93009 


.89795 


.86754 


0. 70 


.09811 


.99623 


.98138 


.96331 


.92872 


. 89609 


. 86529 


0. 75 


.99808 


. 99616 


. 98106 


. 96269 


. 92760 


. 89458 


. 36349 


0. 80 


.99805 


.99611 


.98079 


.96218 


. 92669 


. 89338 


.86207 


0.35 


.99303 


. 99606 


.98056 


.96176 


. 92596 


. 89242 


. 86095 


0. <90 


.99801 


. 99602 


. 98038 


.96142 


.92537 


.89165 


.86007 


0.95 


.99799 


. 99599 


.98023 


.96114 


.92490 


. 89104 


. 85938 


1.00 


.99798 


. 99596 


.98011 


.06092 


. 92452 


. 89056 


.85885 


Z1 


.99856 


.99712 


.98573 


. 97182 


.94501 


.91950 


.89521 


7.BAR 


.99792 


.99585 


.97959 


. 96000 


.92308 


.83889 


.85714 
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ALI 


PHA= 3.0 


P7ER0 


1 = 0.60 
LAMBDA 


MU= 


3.11 


TIME 


.005 


.010 


.050 


.100 


.200 


.300 


0.05 


. 99975 


.99950 


.99750 


.99501 


.99005 


.98511 


C.IO 


.99950 


.99900 


.99502 


.99006 


.98022 


.97048 


0.15 


.9Q925 


.99851 


.99255 


.98516 


. 97054 


. 95614 


0. 20 


.99901 


.99802 


.99012 


.98034 


.96106 


.94217 


0.25 


.09877 


.99753 


.98773 


.97562 


.95184 


. 92865 


0. 30 


.99853 


. 99706 


. 98541 


. 97103 


.94291 


.91563 


0.35 


.99830 


.99661 


.98315 


.96659 


.93434 


. 90320 


C.40 


.99808 


.99616 


.98098 


.96233 


.92615 


. 8Q138 


0.45 


.99787 


. 99574 


.97889 


.95826 


.91837 


.88023 


o 

• 

o 


.99766 


.99534 


. 97691 


.95439 


.91102 


. 86978 


0. 55 


.99747 


. 99495 


.97503 


.95074 


.90413 


. 86002 


0.60 


.99729 


. 99459 


.97325 


.94730 


.89768 


.85097 


0.65 


.90712 


.99424 


.97158 


. 94408 


. 89170 


. 84262 


0.70 


.99695 


. 99392 


.97002 


.94108 


. 88616 


.83495 


0.75 


.99680 


.99362 


.96857 


.93830 


. 88106 


. 82795 


0.80 


.99666 


. 99334 


.96722 


.93572 


. 87638 


.82159 


0.85 


.99653 


.99308 


.96597 


.93334 


.87210 


.81583 


0.90 


.99641 


.99284 


.96481 


.93116 


. 86822 


. 81064 


0.<55 


.99630 


. 99261 


.96375 


.92917 


. 86469 


.80599 


1.00 


.99619 


.99241 


.96278 


.92735 


. 86151 


.80183 


n 


.99782 


.99564 


.97849 


.95767 


.91797 


. 88071 


ZB4R 


.99520 


.99045 


.95399 


.91202 


. 83827 


.77556 



2H 



.400 

.9S020 
. 96083 
. 94195 
.92365 
.90603 
. 88916 
. 87313 
, 85799 
. 8437Q 
. 83055 
.81828 
. 80698 
.79664 
.78721 
. 77868 
.77099 
. 76410 
.75796 
.75251 
.74770 

.84572 

. 72158 



ALPHA= ^.0 P7FR0= 0.70 VU= 1.62 



LA»^R0A 



TI ME 


.005 


.010 


.050 


.100 


. 200 


. 300 


. 400 


0.05 


. 90975 


.99950 


.99750 


. 99501 


.99005 


.98512 


. 98020 


n.io 


. 99950 


. 99O00 


.09502 


.99007 


.98023 


.97049 


.96085 


0.15 


.99925 


,99851 


.99256 


.98518 


. 07059 


. 95621 


. 04205 


O 

• 

o 


.90901 


,99802 


. 99015 


.98040 


.96120 


. 94237 


. 92392 


0.25 


, 9Q077 


. 99755 


.98780 


.97576 


.95212 


.92907 


. 9065R 


0.10 


,90854 


. 99709 


, 98554 


, 97129 


. 94343 


. 91640 


. 8O016 


0.35 


.09832 


.99665 


.98336 


.96702 


.93517 


.90443 


. 97474 


C.40 


.99811 


.99623 


.98130 


.96296 


. 92739 


.80321 


. 860^8 


0.45 


,00791 


.09583 


. 97934 


. 95915 


.92011 


. 88279 


,8471? 


0. 50 


.99773 


, 99546 


.07751 


.95559 


.91335 


.87318 


. 83AQR 


0.55 


,99755 


,99511 


.97581 


.95227 


.90711 


. 3643« 


. 82304 


0.60 


.09739 


.99478 


. 97422 


.94921 


.00140 


.85637 


.81397 


0.65 


,09724 


. 994^8 


,97277 


,04640 


. 89619 


.84913 


. 80503 


0.70 


.09710 


,99420 


.97143 


,94384 


. 80147 


. 84263 


. 70708 


0. 75 


,90697 


. 99395 


. 97020 


. 94150 


.88721 


,83682 


.79004 


o 

or 

9 

o 


.90685 


.99372 


.96909 


.93939 


. 38340 


. 83167 


. 78386 


0,85 


.99675 


,09351 


. 96808 


. 93 748 


, 87909 


. «271 1 


.77R46 


0.90 


.99665 


, 99332 


,96717 


.93577 


.87696 


.82311 


. 77^77 


0,95 


.99656 


.99314 


. 96635 


.93423 


. 87429 


. 81962 


. 76973 


1.00 


. 99649 


. 99299 


.06562 


.93287 


.3719^ 


. 81658 


.76627 


71 


.09791 


. 99583 


.97940 


.959^5 


.92138 


. 38562 


.95200 


7BAR 


.99587 


,99178 


.96021 


.92347 


,85782 


. 80088 


. 75104 



25 



ALPHA= 3.0 PZERD 



TIME 


.005 


.010 


.050 


0.05 


.99975 


.09950 


.99750 


0.10 


.99950 


.99900 


.99503 


0.1 5 


.99926 


.99851 


. 99259 


0,20 


.99902 


.99803 


.99021 


0.25 


.99879 


.99757 


. 98792 


0. 30 


.99857 


.99713 


.08575 


0.35 


.99836 


.99672 


. ?8370 


0.^0 


,99816 


. 99633 


. 98179 


0.45 


.99798 


.99597 


.98002 


0. 50 


.09782 


.99564 


. 97841 


0.55 


.99766 


.09534 


. 97694 


0.60 


.99753 


9506 


.97561 


0.65 


.99740 


. 99482 


.97441 


0.70 


.99729 


. 99460 


.97335 


0.75 


.99719 


.99440 


.97240 


0.80 


.99711 


. 99422 


.97157 


0.85 


.99703 


.99407 


.97083 


0.90 


.99696 


. 99393 


. 97019 


0.95 


. 99690 


.99381 


.96963 


1.00 


.99685 


.99371 


.96914 


Zl 


.99803 


. 99608 


. 98062 


ZRAR 


.99652 


. 99307 


. 96629 



0.80 


MU= 


4. 30 




LAMBDA 


.100 


.200 


.300 


.400 


99501 


.99005 


.98512 


.98021 


99008 


.98025 


. 97052 


. 96089 


98523 


.97067 


.95634 


.94221 


98052 


.P6142 


. 94270 


.92436 


97600 


.95259 


.92976 


.90750 


97171 


.94426 


.91762 


.89177 


96768 


.93648 


. 90636 


. 87723 


96394 


.92930 


.89603 


.86407 


96049 


.<92274 


.88665 


. 85215 


95735 


.91678 


. 87820 


.84151 


95450 


.91143 


. 87068 


.83210 


95193 


.90666 


. 86402 


. 82385 


94964 


. 90243 


. 85818 


. 81668 


94760 


.89872 


.85309 


.81050 


94581 


. 89547 


. 84869 


.80522 


94423 


.89264 


.84491 


.80074 


94284 


.89020 


. 84169 


. 79697 


94164 


.88810 


. 83396 


.79383 


94060 


.88631 


.83666 


.79122 


93970 


. 88479 


. 83474 


. 78908 


96183 


. 92595 


. 89218 


. 86038 


93478 


. 87755 


. 82692 


.78182 



26 





ALi 


5HA= 3.0 


PZERO 


= 0.90 


MU= 


5.35 












LAMBDA 








TIMP 


.005 


.010 


.050 


. 100 


. 200 


. 300 


.400 


0.05 


.99975 


, QQ05O 


. 99750 


.99502 


.99006 


.98512 


. 98021 


0.10 


.90950 


. 99901 


.99504 


,99010 


. 98029 


.97058 


.96097 


0.15 


.99926 


.90852 


.99263 


,98531 


. 97084 


. 95659 


.94255 


0.20 


. 90903 


. 99806 


. 90032 


.98074 


.96187 


.94337 


.92524 


0.25 


.99381 


.99762 


.98815 


.97646 


. 95349 


.93110 


. 90927 


0. 30 


.99861 


. 99721 


. 98614 


. 97250 


. 94581 


. 91992 


. 89479 


0. 35 


.99842 


. 99684 


.98431 


.96890 


.93887 


.90988 


. 88189 


0.40 


.99825 


.99650 


.98266 


.96566 


. 93268 


. 90100 


. 87055 


0.45 


.99510 


. 99620 


.98119 


. 96280 


.92724 


.89324 


.86074 


0.50 


.99797 


.99594 


.97990 


.96029 


.92251 


. 88656 


. 85236 


0. 55 


.99785 


. 99571 


. 97877 


. 95810 


.91843 


. 88087 


. 84528 


0.50 


.99775 


.99551 


.97779 


.95623 


.91496 


. 87606 


. 83938 


0.65 


.9Q766 


. 99533 


.97695 


.95462 


.91203 


. 87205 


. 83450 


0.70 


.99759 


. 99518 


,97623 


.95326 


.90957 


.86873 


.83052 


0.75 


.99752 


.99506 


.97563 


.95211 


.90753 


. 86601 


. 82731 


C.80 


.99747 


.99495 


.9751 2 


.95115 


. 90585 


. 86379 


.82474 


0. 55 


.Q9742 


. 99486 


. 97469 


. 95036 


. 90447 


. 86201 


.82971 


0. 00 


.99939 


. 99478 


.97433 


.94970 


.90334 


. 86059 


.32112 


0.95 


.99735 


. 99472 


. 97403 


. 94915 


. 90244 


. 85047 


.81990 


1. 00 


.99733 


. 99466 


,97379 


.94871 


.90171 


.85858 


.81896 


Zl 


.99822 


. 996 44 


.98241 


.96534 


.93266 


. 90181 


. 87268 


ZBAR 


.99720 


.99442 


.97273 


, 94690 


. 80916 


. 85600 


. 81679 



27 



ALPHA= -3.0 



PZER0= 0.95 



MU= 6.30 



LAMBDA 



TIME 


.005 


,010 


.050 


. 100 


. 200 


. 300 


.400 


0.05 


.99975 


.99950 


.99751 


.99502 


.99006 


,98513 


.98022 


0. 10 


.99950 


.99901 


.99505 


.99012 


.98034 


.97066 


, 96107 


0.15 


.99927 


. 99853 


.99268 


.98541 


.97105 


,95689 


.94295 


0.20 


.^9904 


.99808 


.99045 


.08100 


. 96238 


, 94413 


.92624 


0.25 


.99883 


. 99767 


. 98841 


.97696 


.95448 


.93257 


.91119 


0.30 


.99865 


.99730 


.98656 


.97333 


,94744 


.92233 


. 89796 


0.35 


,99848 


. 99697 


. 98493 


. 97013 


. 94129 


,91345 


.88656 


0.40 


. 99834 


.99668 


.98352 


.96736 


,93600 


.90537 


.87692 


0.^5 


.99821 


,99643 


.98230 


, 96500 


.93153 


. 89952 


. 86891 


C.50 


.99811 


. 9962? 


.98127 


,96300 


.92779 


.89426 


,86235 


0.55 


.99802 


.99604 


.98041 


.96134 


,92470 


, 88997 


. 85705 


0. 60 


.99794 


. 99589 


.97969 


.95997 


,92219 


.88652 


,85284 


0.65 


.99788 


.99577 


.97910 


.95884 


.92015 


,88377 


. 84954 


0.70 


,99783 


. 99567 


,97862 


.95793 


,91853 


. 88160 


. 84698 


0.75 


.99779 


. 99559 


.97823 


.95720 


.91724 


,87991 


.84502 


0.80 


.99776 


,99552 


.97792 


,95661 


.91623 


.87861 


. 84354 


0.85 


,99773 


.99547 


,97766 


.95615 


.91544 


. 87762 


.84245 


0.90 


.09771 


. 99542 


.97746 


.95578 


.91484 


.87688 


.84164 


0, 95 


.09769 


.09539 


.97730 


,95549 


.91437 


. 87632 


. 84106 


1.00 


.99768 


.99536 


.97718 


.05527 


.91402 


.87591 


.84065 


Z1 


.99837 


. 99674 


. 98388 


.96821 


.93815 


.90970 


.88274 


ZBAR 


.99762 


,99526 


.97674 


.95455 


.91304 


.87500 


.84000 



28 



400 





ALI 


o 

• 

II 

< 

X 


PZERO 


= 0.60 


MU= 


4. 20 












LAMBDA 








T IMF 


.005 


.010 


.050 


. 100 


.200 


.300 


.400 


C. 05 


. 00975 


. 99950 


.99750 


.99501 


.99005 


.98511 


.98020 


C.IO 


.99950 


, QQ9on 


.99501 


.09005 


. 98020 


.97045 


, 96080 


0. 1 5 


.00925 


.99850 


.99253 


.98512 


.97047 


.95604 


.94182 


0. 20 


.99900 


.99801 


.99007 


.98025 


. 96088 


.94191 


.92330 


0.?5 


.99876 


.00752 


.98764 


.97544 


, 95148 


.92812 


.90533 


0. 5^ 


. 90^5? 


. 99703 


.98525 


.97072 


.04231 


.91474 


. 8870Q 


C.25 


.90828 


.09656 


.93292 


.96613 


. 93343 


.90185 


. 87137 


0.40 


.09805 


.99610 


.98065 


. 96168 


, 92488 


.88052 


, 85556 


0. 45 


.09782 


.99565 


.97847 


.05741 


.91671 


.87781 


. 84064 


0,50 


.99761 


.09523 


.97637 


.95334 


. 90897 


. 86678 


. 82667 


0. 55 


.99741 


. 99482 


. 97430 


.94948 


.90168 


.85647 


.81371 


0.60 


.09721 


.99444 


.97251 


.94584 


. 89488 


. 84601 


. 80177 


0.65 


.00703 


. 99407 


.97075 


.94245 


. 88856 


. 83811 


. 79087 


0. 70 


.99686 


. 99373 


.96911 


.93929 


.88274 


.83007 


.78101 


0.75 


.09670 


.99342 


.96758 


.93638 


.87742 


. 82278 


.7721 4 


0.80 


. 99656 


.99313 


.9661 8 


.93370 


.87258 


.81622 


. 76425 


0.85 


.90642 


.09286 


. 96489 


.93126 


. 86821 


.81036 


. 75727 


0.00 


.00630 


.99261 


.96371 


.02905 


. 86428 


. 80515 


.75117 


0. 9 5 


.99618 


. 99238 


.96264 


.92 704 


. 86078 


. 80057 


.74586 


1.00 


.90608 


.9921 8 


.96168 


.92524 


.85767 


. 79657 


. 74130 


71 


.09776 


.99553 


.97795 


.95661 


.91597 


. 87736 


. 84212 


ZRAR 


.99526 


. 99057 


.95455 


.91304 


.84000 


. 77778 


. 72414 



29 



ALPHA= 4.0 



PZER0= 0.70 
LAMBDA 



MU= 4.76 



TIME 


.005 


.010 


.050 


. 100 


.200 


. 300 


.400 


0.05 


.99975 


.99950 


. 99750 


.99501 


.99005 


.98511 


.98020 


0.10 


, 99950 


.99900 


.99501 


.99005 


.98020 


.97045 


. 96080 


0.15 


.99925 


.99850 


.99254 


.98513 


. 97049 


. 95606 


. 94184 


0.20 


. 99900 


. 99801 


.99009 


.98027 


.96093 


.94198 


.92340 


0.25 


.99876 


.99752 


.98767 


.97550 


.95161 


.92831 


. 90553 


0. 30 


.99852 


.99705 


.98532 


. 97086 


.94258 


.91513 


.88851 


0. 35 


.99829 


.99658 


.98304 


.96637 


.93390 


.90256 


. 87229 


0.40 


.99807 


.99614 


.98085 


.96208 


.92565 


. 89066 


. 85705 


0.45 


.99785 


. 99571 


.97876 


. 95800 


.91787 


.87952 


.84287 


0.50 


.99765 


.99531 


.97680 


.05418 


.91061 


. 86919 


.82982 


0.55 


.99746 


.99494 


.97496 


. 95061 


.90389 


. 85971 


.81791 


0.60 


. 99729 


. 99459 


.97325 


.94731 


.89773 


. 85108 


. 80718 


0.65 


.99713 


.99426 


.97168 


.94429 


.89213 


. 84330 


. 79759 


0.70 


.99698 


. 99396 


.97025 


.94153 


. 88708 


.83636 


.78911 


0.75 


.99684 


. 99369 


.96894 


.93904 


.88256 


. 83021 


.78169 


0.80 


.99672 


.99345 


.96776 


.93681 


. 37855 


. 82481 


.77526 


0.85 


.99661 


. 99323 


.96671 


.93482 


.87501 


.82012 


.76974 


0.90 


.99651 


.99303 


.96576 


.93305 


. 87192 


. 81608 


-76506 


0.95 


.99642 


. 99285 


.96493 


.93150 


. 86924 


.81263 


.76113 


1.00 


.99634 


. 99269 


.96419 


.93013 


.86692 


.80970 

# 


. 75787 


Z1 


.99784 


. 99568 


.97870 


.95807 


.91377 


.88191 


. 84730 


Z8AR 


.99582 


.99167 


.95968 


.92248 


.85612 


. 79866 


. 74843 



30 



ALPHA= 4,0 P/FRO 



TIME 


.005 


.010 


.059 


0.05 


.9^975 


,09050 


.99750 


0. 10 


. 09950 


. 09900 


.09502 


0.1 5 


.90075 


. 09950 


.99755 


0. ?0 


.99901 


.99801 


.90011 


0. 25 


.09877 


• 0975^ 


.98773 


0.30 


.09853 


.99707 


. 00544 


0. 3 5 


.99831 


. 0966? 


.98324 


O.AQ 


.99810 


.09620 


.98118 


0.45 


.99790 


. 90581 


,07925 


0. 50 


.99772 


. 99545 


.97747 


0.55 


.99755 


.90512 


,07585 


0.60 


.99740 


. 00481 


.97438 


0.65 
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APPENDIX B 



EXAMPLE PROBLEM 

Suppose a system consisting of 15 components is connected 
in series. Each component , operating independently of the 
others, has an exponential failure density with parameter X 
and a special Erlang repair density with parameters y and 
a. Assume that uptimes and downtimes are all independent. 

a. What is the average availability of the system over 
the first mission time if; 

(1) X = .005 failures per mission time, y = 4.8 
repairs per mission time, and a = 2? 

(2) X = .40 failures per mission time, y = 4.2 
repairs per mission time, and a = 4? 

b. What is the long-run average availability of the 
system if; 

(1) the parameters are the same as in a(l)? 

(2) the parameters are the same as in a(2)? 

Solution : 

a. 

(1) The average availability of a single com- 
ponent over the first mission time is Z1 = .99856 (from 

Appendix A). Letting n be the average system availability 

s 

over the first mission time; 

rig « (Zl )^5 = (.99856)^5 = . 9786 . 
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(2) In this case Zl=. 84212 and n =(.84212)^^= .0760 . 

s 

b. 

(1) The long-run average availability of a single 
component Is ZBAR = .99792. 

Letting n be the long-run average availability of the 
system, 

W = (ZBAR)^^ = (.99792)^5 = .9693 • 

This compares with the average system availability over 
the first mission time of .9786. 

(2) In this case ZBAR = .72414 and n". = .0079 . 

This compares with the average system availability over 

the first mission time of .076. 

The above comparisons Indicate the use of the long-run 
average availability to approximate the system availability 
over the first mission time Is a conservative practice. 

These cases and others also Indicate that the lower the 
value of ZBAR the more conservative Is the practice. 
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APPENDIX C 



FORTRAN IV COMPUTER PROGRAM FOR THE 
INVERSION OF n*(s) BY NUMERICAL METHODS 



THIS PROGRAM COMPUTES VALUES OF SYSTEM AVAILABILITY BY 
INVERTING THE LAPLACE TRANSFORM OF THE AVAILABILITY 
FUNCTION, 

IMPLICIT REAL*B( A-H),R5AL*fl(0-Z > 

C0MPLEX*16 S» ETA, EX PON, CONST, TOP ♦bottom 
DIMENSION XCnF(?OI ,COF(20) ,R0nxR(20) , ROOT I ( 20 ) , S( 20 ) , 
1 FXP0N(20 ) , RETA(20,7 ) ,YLAMDA( 7 ) , ZBAR ( 7 ) , Z I ( 7) ,ETA( 20) 
DATA YLAMOA/O.Or'5,0.010,0.05,0.1 ,0 . 2 , 0 . 3 , 0 .4/ 

JCASF=0 

ICASE IS THF NUMBER OF CASES TO RE PROCESSED THIS RUN, 
PEAD(5, 1 ) ICASE 
1 FORMAT(I2) 

E REAO(5, 10)ALPHA,XMU, PZERO 
10 FORMAT(F3,0,2F10,3) 

JCASE = JCASF4-1 

WRITE(6,6) JCASE 

6 FORMAT! IHI ,45X, 'CASE* , 14) 

DO 800 K=l,7 
XLAMDA=YLAM0A(K) 

WRITE(6,7) 

7 F0RMAT(//10X, • INPUT DATA*) 

WRITE(6,8) ALPHA, XMU,XLAMDA 

8 FORMAT ( /lOX, • ALPHA = *,F3,0,5X,* MU =*,F10,3,5X, 
I'LAMBOA =',F10,3) 

PUT ZEROS IN ALL ARRAYS, 

DO 80 L=l, 20 
XC0F(L)=0,0 
C0F(L)=0,0 
ROOTR(L)=0.0 
ROOTKL )=0,0 
S (L) = (0, 0,0,0) 

FTA(L)=(0, 0,0,0) 

80 EXPONCL )=( 0,0,0,0) 

C 

M=ALPHA 
N = M+1 

IF(M,NE, ALPHA )G0 TO 40 
C 

C COMPUTE COEFFICIENTS OF THE POLYNOMIAL, 

CALL COEFF ( AL PH A, X MU, X L AMO A , XCO F ) 

C 

C find the ROOTS OF THE POLYNOMIAL, 

C DPOLRT IS A PROGRAM LIBRARY SUBROUTINE. 

CALL DPOLRT(XCOF,COF,N,ROOTR,RnOTI , lER ) 

C 

60 WRITE(6,61 ) lEP 

C lER IS AN ERROR CODE FROM SUBROUTINE OEPOLRT, 

61 FORMAT(//10X, *IER=*,I2) 

WRITE! 6,18) 

IB FORMAT! /13X, ‘ROOTS OF THE POLYNOMIAL,*) 

WRITE(6 19) 

19 FORMAT! /15X, 'REAL* ,6X, • IMAGINARY'/ ) 

14 WRITE!6,16) !ROOTR! I ),ROOTI! I), I=l,N) 

16 F0RMAT!10X,2F12,5) 

C 

C CHECK FOR MULTIPLE ROHTS, 

CALL CKMULT !R00TR, ROOT I,M) 
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c 

C PUT THE ROOTS INTO THE COMPLEX ARRAY S(20). 

00 100 1=1, N 

100 Sm =OCMPLX(ROOTR( I ) ,R00TI ( I ) ) 

BEGIN INVERSION BY HEAVISIDE'S EXPANSION FORMULA, 

00 150 1=1, N 
T0P=XMU+S( I ) 

RnTT0M=Tnp + M* (XL AMOA+S ( I ) I 
C0NST=T0P/R0TT0M 
DO 150 J=l,20 

FXP0N( J )=C0NST*C0EXP( J*S( I)*0. 05> 

150 ETA( J)=ETA( J) +eXP0N( J ) 

WRITE! 6,160) 

160 FORMAT (///lOX, 'VALUES .OF THE AVAILABILITY FUNCTION AT' 
I' .05 MISSION-TIME INCREMENTS.'/) 

DO 170 1=1,20 
T = 0.05* I 

170 WRITE(6,171 )T,FTA( I ) 

171 FORMAT! 10X,F4. 2,2F10.6) 

SUM=0. 0 
00 155 L=l,20 
RETAIL, K) = PEAL(FTA!L) ) 

155 SUM=SUM+RETA! L ,K) 

TOTSUM=SUM+SUM+l,-RETA! 20,K) 

71 IS THE AVERAGE AVAILABILITY OVER THE FIRST MISSION-TIM 
Z1!K)=TOTSUM/40. 

ZRAR IS THE LONG RUN TIME AVERAGE AVAILABILITY. 

900 ZBAR(K )=XMU/( ALPHA*XLA MOA+XMU) 

CALL OUTPUT!YLAMOA, ALPHA,PZERG,XMU, Zl, 7BAR ,RETA) 

30 IF! JCASE-ICASE)5,50,‘?0 

40 WRITF(6,41) 

41 FORMAT! lOX , 'ERROR- ALPHA MUST RE AN INTEGER') 

GO TO 30 

50 STOP 
END 



SUBROUTINE CK MULT ! ROOTR , ROOT I , M ) 

IMPLICIT REALMS! A-H) ,REAL*8( 0-Z) 

DIMENSION S (1 ), ROOTR! 1 ) ,ROOTI! 1 ) 

THIS SUBROUTINE CHECKS FOR MULTIPLICITIES AMONG THE ROOTS 

MULT=1 
DO 50 1=1, M 
K=M+1-I 
DO 50 J=1,K 

50 IF! ROOTR! I ) .FQ. ROOTR ( I+J ). AND. ROOT I ! I ) ,E0. ROOT I ( I+J) ) 
1MULT=MULT+1 

IF!MULT.NE. 1 ) GO TO 60 
C 

WRITE!6,51 ) 

51 FORMAT! ///lOX, ' THERE ARE NO MULTIPLE ROOTS.') 

RETURN 

C 

60 WRITE(6,61) 

61 FORMAT !///10X, 'THE POLYNOMIAL HAS MULTIPLE ROOTS, THIS 
I'PROGRAM WILL NOT ACCURATELY HANDLE THIS CASE,') 

RETURN 

END 
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SUBROUTINE COEFF ( A L PHA , XMU , XLAMnA,XCQF) 

IMPLICIT REAL’«'8( A-H),RFAL^8(0-Z I 
DIMENSION XCOFd) 

THIS SUBROUTINE COMPUTES THE COEFFICIENTS OF THE EXPANDED 
POLYNOMIAL IN THE DENOMINATOR OF THE LAPLACE TRANSFORM. 

M=ALPHA 

N2=M+2 

X=ALPHA+1.0 

XCOFd )=0.0 

XC0F(M4-2)=1.0 

CHOOZ=M 

DO 20 1=1, M 

BKT=XLAMDA+( XMU*n /( X-II 
XCOFI I+l ) = CHOOZ + BKT*XMU*»(M-n 
20 CHooz=CHOoz-*=(M-n/ ( di ) 

THE COEFFICIENTS ARE NOW IN THE ARRAY, XCOF, 

WRITE(6,17) 

17 FORMAT ( //lOX, 'COEFFICIENTS OF THE POLYNOMIAL') 

WR ITEI6, 12) (XCOFd ) ,1=1, N2 ) 

12 FORMAT! /lOX, 10F12.6) 

RETURN 

END 



SUBROUTINE OUTPUT ( YLAMDA , ALPHA , PZE RO , XMU, Z 1 , ZBAR , RETA ) 
IMPLICIT REAL’''0(A-H),REAL*8(O-Z) 

DIMENSION YLAMOA( 1 ) , Zl( 1 ) , ZB AR ( 1 ) , R ETA ( 20 , 7 ) 

WRITEC6 ,1 ) ALPHA, PZERO, XMU 

1 FORMAT! 1H1,////////25X,'ALPHA=',FA,1,5X,'PZER0=«,F5.2, 
15X, ' MU= ' ,F5.2 ) 

WRITE!6,11 ) 

11 FORMAT! //46X, ' LAMBDA' ) 

WRITE (6, 2) !YLAMOA!L ),L=1,7) 

2 FORMAT! /22X,FA. 3,6!AX,F4.3 ) ) 

WRITE!6,12 ) 

12 F0RMAT!15X, 'TIME'/ ) 

DO 50 1=1,20 
T=0.05*I 

50 WRITE(6,3)T,!RETA!1,J),J=1,7) 

3 FORMAT! 15X,F4.2,7!2X,F6. 5) /) 

WRITE!6,4) !Z1( I), 1=1, 7) 

4 FORMAT !/ 16 X,'Z1' ,1X,7!2X,F6.5)//) 

WRITE!6,5) ! ZBAR! I) , 1=1,7) 

5 F0RMAT!15X , ' ZBAR' ,7!2X,F6. 5 ) ) 

RETURN 

END 
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